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o\ ■ 

' We study finite size effects for the gap of the quasiparticle excitation spectrum in the weakly 

, interacting regime one-dimensional Hubbard model with on-site attraction. Two type of corrections 

' to the result of the thermodynamic limit are obtained. Aside from a power law (conformal) correction 

^ ' due to gapless excitations which behaves as 1/Na, where Na is the number of lattice sites, we obtain 

^ , corrections related to the existence of gapped excitations. First of all, there is an exponential 

■ correction which in the weakly interacting regime {\U\ <C t) behaves as ~ exp(— A^aAoo/4t) in the 
' extreme limit of Na^ac/t ^ 1, where t is the hopping amplitude, U is the on-site energy, and Aoo is 

00 \ the gap in the thermodynamic limit. Second, in a finite size system a spin-fiip producing unpaired 

fermions leads to the appearance of solitons with non-zero momenta, which provides an extra (non- 
exponential) contribution 5. For moderate but still large values of NaAoa/t, these corrections 
D ' significantly increase and may become comparable with the 1/Na conformal correction. Moreover, 

in the case of weak interactions where Aoo ^ t, the exponential correction exceeds higher order 
power law corrections in a wide range of parameters, namely for Na < (8t/ A^a) ln(4t/|f7|), and so 
does 5 even in a wider range of Na. For sufficiently small number of particles, which can be of the 
, order of thousands in the weakly interacting regime, the gap is fully dominated by finite size effects. 

^ , I. INTRODUCTION 

a , 

O . The one-dimensional (ID) Hubbard model with only on-site interaction is exactly solvable by the Bethe Ansatz, 
and the properties of this model have been widely studied Low lying excitations of the Hubbard model with 
attraction are gapless charge excitations and gapped spin excitations, whereas for the repulsive model the gap exists 
CO ' only at half filling in the charge sector. 

^ [ Finite size corrections to the ground state energy Eq, due to the gapless part of the spectrum, follow from conformal 
field theory 0, H, 13, Hi and have the form: Eq — eoL = —TTvh/dL, where L = NaQ is the size of the system, Na is 
the number of lattice sites, a is the lattice constant, eo is the energy per unit length in the thermodynamic limit, 
and V is the velocity of gapless excitations. Finite size corrections to the energies of low-lying gapless excitations are 
also proportional to v /Na, and the proportionality coefficient depends on the scaling dimensions of the primary fields. 
However, finite size corrections originating from the gapped sector remained unknown for the Hubbard model. Finite 

■ size effects are expected to be important for sufficiently small systems, such as cold atoms in a ID optical lattice, 
\ where the number of particles and lattice sites ranges from several tens to several hundreds 0] ■ 

IL^ ■ In this paper we solve the Bethe Ansatz equations for a finite number of particles and calculate finite size corrections 
to the gap for the attractive Hubbard model. As expected, there are power law corrections due to gapless 

excitations, and we also find contributions related to the existence of gapped excitations. First of all, there is 
?H ' an exponential correction which in the weakly interacting regime {\U\ ^ t) behaves as ~ exp(— iVaAoo/4t) in the 
extreme limit of NaA^c/t ^ 1, where t is the hopping amplitude, U is the on-site interaction, and A^o is the gap in the 
thermodynamic limit. Second, in a finite size system a spin-flip producing unpaired fermions leads to the appearance 
of solitons with non-zero momenta, which provides an extra (non-exponential) contribution. For moderate but still 
large values of NaAoo/t, these corrections may become comparable with the 1/Na conformal correction. Moreover, in 
the case of weak interactions where Aoo ^ t, the exponential correction exceeds higher order power law corrections 
in a wide range of parameters, namely for Na < (8t/Aoo) ln(4t/|J7|), and so does S even in a wider range of Na. We 
find that the value of the gap increases with decreasing the system size and show how the gap becomes dominated by 
finite size effects in sufficiently small systems. 

The paper is organized as follows. In section II we introduce the Hubbard model together with related Bethe Ansatz 
equations, and discuss the thermodynamic limit. In section HI we present a general approach for finding finite size 
corrections to the ground state energy and to the gap and discuss the structure of the gap. Section IV contains our 
results for corrections due to the gapped sector at half filling, and section V the results for power law corrections. In 
Section VI we discuss our numerical and analytical results, and in Section VII present for completeness a perturbative 
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approach for solving the Bethe Ansatz equations in the limiting case of L <C at/U . In Section VIII we conclude. 



II. GENERAL EQUATIONS, THERMODYNAMIC LIMIT 

The Hubbard model for a system of interacting spin-1/2 fermions on a lattice is described by the Hamiltonian 

H = -t (4<TCi+i,a + 4+i.<TCj:'^) + ^$I"j.T'^:i-,i. (1) 

'T=Ta;j=i j=i 

where the subscript j labels the lattice sites. The index a labels the spin projection, cj ^ and Cj^a are the creation and 

annihilation fermion operators, and rij^a- — c^j a'^j.a are the particle number operators. Below we express all quantities 
having the dimension of energy in units of t, and quantities having the dimension of length in units of the lattice 
constant a. 

Lieb and Wu have solved the Fermi-Hubbard model by means of the Bethe Ansatz . The corresponding eigenvalue 
equations read 



^2arctan — smkj\^ ^ ^^j^ ^_ ^ 2 arctan (^-^y— > a = 1, ...M, 



(2) 



M 



Nakj = 2nl, ~ Yl 2 arctan [ ^'^ , J = 1, -N, (3) 

where u = |C/|/4t, M is the number of spin-down fermions, and N is the total number of particles. The energy of a 
given state is expressed through the charge momenta kj-. 

N 

En = -2^(^08 kj, (4) 

and it depends on the spin rapidities Aq only implicitly through the coupled equations ([2]) and ((S)). For the ground 
state the quantum numbers Ij and Ja are integers or half-odd integers depending on the parities of N and M: 



N + M + 1, _ , ^ M , 

T 



Ja = ^ (mod 1), Ij = —(mod 1). (5) 



For the Hubbard model with attraction (U < 0), there is a gap in the spectrum of spin excitations. Considering a 
finite number of particles we define the gap A as: 

2 A = £;Ar+2(^T + 2, iVi, [/) + En-2{N^ ~2,N^,U)- 2En{N^, N^,U), (6) 

where i?Ar(iV|, iVj^; U) is the ground state energy for a system with spin-up and Ni = N — N]- spin-down fermions 
in a lattice with Na lattice sites, at the interaction strength U. This definition is convenient as it does not change the 
parity of the quantum numbers Ij and Ja- Without loss of generality, we may put < N-^. 

In the thermodynamic limit, where N 00 and Na 00 while keeping constant densities n — N/L and = N [/L, 
Eq. (O leads to the same result as the definition A En^ (iV| + 1 , iV| ; C/) - 2En^ {N^ ,Ni;U) + Em^ (7V| - 1 , iV| ; [/) 
introduced in Refs. [j. 

In the thermodynamic limit the density of momenta k and the density of rapidities A are defined as p{k) = L^^dl/dk 
and cr(A) = L~^dJ/dX, respectively. Then, the Bethe Ansatz equations for the ground state of the repulsive model 
become: 

1 2u 1 ['^ u 

aoo(A) + - j ^ _^_^^_^aoo(A')dA' = - j ^ -^-^^-—^p^ik)dk, (8) 

where the constants B and Q are given by 

N f'^ Ni 

— =n=/ p^{k)dk, croo(A)dA, (9) 
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and the ground state energy follows from the relation: 

= -2Na / /9oo(fc)cos/e dk. (10) 

J-Q 

The value of the spin gap for an arbitrary filling factor has been calculated in the thermodynamic limit in Refs. 
[3, ■ In the case of weak attraction the result is 

_ 16sin^/^(7rn/2)VM / tt sin (7rn/2) \ 

^oc — GXp I - I , (11) 

TT \ 2u J 

and the validity of Eq. (jlip requires large values of the exponent. 

The particle-hole symmetry and the symmetry with respect to interchanging spin-up and spin-down fermions allow 
one to establish relations between ground state energies of the repulsive and attractive Hubbard models Q: 

E{N^,N^- U) = -{Na - A^T - N!^)U + E{Na -N^,Na~N^-U) 

= N^-U + E{N^,Na - Ni;-U) = N^U + E{Na ~ N^,Ni;~U). (12) 

Using Eqs. ^ and (|12p we can express the spin gap for the attractive Hubbard model through the energies of the 
repulsive model: 

2A = 2\U\ + En^-2{M ~2,Na- M, \U\) + EN^-2{Na -M-2, M, \U\) - 2EnAM, Na - M, \U\). (13) 
For the half- filled case {N = 2M ^ Na), with N^ = N^ = M, Eq. ^ takes the form: 

A=\U\+ EN^^2{Na/2, Na/2 - 2, \U\) - ENANa/2, Na/2, \U\). (14) 

Below we calculate the gap for the attractive Hubbard model. For this purpose we first perform calculations of 
ground state energies for the repulsive Hubbard model, where the momenta kj are real numbers, and then obtain A 
for the attractive model by using Eq. (|13p (Eq. (14) for the half filled case). This allows us to find exponential finite 
size corrections, which is not possible in direct calculations for the case of attraction where kj are complex and can 
be found only with an exponential accuracy. 



III. FINITE SIZE CORRECTIONS. GENERAL APPROACH 

Thus, in order to calculate finite size corrections to the gap we have to obtain the three energies of a finite size 
system, entering the right hand side of Eq. (fT3|) . We will follow the scheme proposed by de Vega and Woynarovich 
[l(j | , which introduces a formalism allowing us to use the Bethe Ansatz in order to calculate finite size corrections to 
the energy of the ground state. The scheme consists of writing the Bethe Ansatz equations ([2|) and ([3]) in the form: 

Z'{X) = ±y iarctan^^^^-— y iarctan^^, Z^(A„) = ^, (15) 
Z^(fc) ^ A+ 1 f;iarctan^^^^^^, Z^ik,) ^ ^. (16) 

2Tr Na TT U ' Na 

Oi 

We then define the densities of momenta k and rapidities A for a finite size system as 

° j=l " /3=1 

dZ'^ 1 1 cos A: ^ X /,„s 

Q — 1 

where Ki{x) ~ 2u/{v? + x^), and K2 = 4ii/(4w^ -f- x^). The densities satisfy the relations: 

crAr(A)aA = ■ 

A_ 



aN{\)d\=—, / p^(k)dk^—, (19) 
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where Q±, A± are determined from the equations 

We first perform calculations for the half-filled case. According to Eq. ([T^ we have to calculate 
En^ {Na/2, Na/2] \U\) and En^ {Na/2, Na/2 - 2; \U\). In the former case we have N ^ Na and N^=Ni^ Na/2 = M, 
and the quantum numbers for the ground state are 

, M - 1 M -\, 

Jo. = { ^,...,-1,0,1,...,^—}, 

With N = Na and M = 7Vq/2, from Eq. ^ we have 



Then, from Eq. we obtain 



J -iW^. / -i!^ (22) 

'-'max — n ' ^max — c\ " V / 



= Z,(Q+) = i (23) 



A+ = c5o; Q+ = TT. (24) 

When calculating the energy En^ {Na/2, Na/2 - 2; \U\) we have N ^ Na-2, = Na/2, and A^j^ = Ara/2 - 2. It is 
convenient to introduce two additional particles with momenta kh = {ki, fcjv„} and spin rapidities = {Ai, Ajv„/2} 
in order to satisfy the conditions Q± — ±7r and A± = ±oo. The discrete Bethe Ansatz equations in this case read 

El Xa — sm kj \ \ ^ I Xa — Aa 

2 arctan I I = 27r + 2 arctan ' 



and Eqs. ^ and HH) lead to 



j=l ^ ^ /3=1 



2m 



^2arctan(^^^^^) -^2arctanf^^i-^), a = 1, ...Ara/2, (25) 



^ u J ^-^ \ 2u 

Na/2 



Nakj = 27r/,- - 2 arctan ( ^!!L^2 ^ ) + 2 arctan ( ^ ) ^ j = 1, ...A^a- (26) 

^1 V " ^ tr V ; 

Note that we added only 4 additional equations for defining the numbers kh, \h- Other equations are exactly the 
Bethe Ansatz equations for Na — 2 particles. The sets J a, Ij are the same as for the ground state of Na particles (|2ip . 
For the case A^ = Na we rewrite equations ((TT)) and p8| in the form of integral equations for the densities (TjVa(A), 

1 r 1 f°° 

o'JV„(A) ^ I Ki{X - sink)pNa{k)dk - — / K2{X - ti)crNa{f^)dfJ. 

1 r 1 f°° 

+ — Ki{X~ sin k)X%^{k)dk-— K2{X- fi}X'j^^ {fi)dfi, (27) 
^TT Zn 

1 1 f°° 

Pn {k) — 1 / ATi (sin fc — A) cos /c {X)dX 

2tt 27r J^^ 

1 /"^ 

+ — cosfc A:i(sinfc- A)X^ (A)dA, (28) 
2tt , _^ 
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with the quantities and X'f^ defined as 
For N — Na — 2 we do the same and obtain 

crAr„-2(A) = — 



1 
1 



Na 
'Na/2 



1 

2^ 



OO 



(29) 
(30) 



Ki{X — sink)pNa~2{k) dk 
+ — Xi(A-sinfc)X]v^_2(fc)dfc-— / K2{X-^i)X'r^^^^i^i)d^l 
- 2^E^i(A-Binfc.) + ^E^^(A-A.), 

1 1 /""^ 

PNa-2{k) = h — / iiri(sinfc - A)crAr^_2(A) cosfc dA 

1 Z""" 1 

— / ifi(sinA: — A)X^ _2(A) coskdX—- — — - 7 iiri(sin A: — A/t) cos fc. 
27r 7_oo " 27rA'a 



The functions X^'^_2 are given by 



(31) 



(32) 



^^„-2(A) = ^ E '5(A - A„) - a^r_2(A), 

^ Q — 1 

. Na 

Xk-2ik) = TT E '^(^ - - /^A'«-2 



(33) 



(34) 



where the summation over kj, Aq, includes the additional numbers fc/i, A/j. 

We first consider the thermodynamic limit where the terms X^^^ in Eqs. (|27p . PT|) and ([5^ vanish. For the 

half-filled case of = Na the solution is known 0] : 



PoC,Na{k) 

Coo, Wo (A) ^ TT f 



1 cos fc Z""" cos(w sin /c) exp(— ww) 



2tt 2tt jo 
coso^A 



cosh WW 



cosh WW 
Jo{ui)dLo. 



Jo(w)(iaj, 



Integrating the density of momenta over dk we obtain 



,N, 



k 1 
(^)-2^ + 2^ 



sin(a; sin fc) exp(— wit) 



uj cosh WW 



Jo(a;)(iw. 



(35) 
(36) 

(37) 



For the case of = iV^ — 2 the solution is easily found 

<^oo,Na.-2 = O-oc.N, 



1 X ^ TT 

}_sr^ f°° cos[w(A - A/0] 



2'KNa ^ 2u cosh[7r(A — sinfc;i)/2M] 



2-KNa — JO 
Ah 



cosht 



exp(— wm) cZw, 



(38) 



PoO,Na-2 = PoO,Na 



COS fc 



COS fc 

27r7V, 



E 



COS [w (sin fc — sin fc/j) 
cosh WW 



exp(— ww) duj 



y- 

^ Oil 



2-KNa ^ 2u cosh[7r(sinfc — Xh)/2u] ' 



(39) 
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and 



2ttN, 



■E 



sin [a;(sin k — sin k^)] exp(— wm) 



2ttN, 



2 arctan j^tanh ^ — (si 



OJ cosh 

sin k — \h) 



duj 



(40) 



For calculating the finite size corrections we should find the differences between the densities PNa.{k), f^NaWi 
<JNa-2{^), PNa-2{k) and their thermodynamic limit values. Subtracting equations of the thermodynamic limit from 
Eqs. ^7}i, (111]) for the case iV = iV^ we obtain: 



^f^N^W = O-atJA) - (Too,iV„ (A) 



2tt 



Ki (A — sin k)SpN^ (k) dk 



2tt 



K2{\ - p)6aN^{p) djjL 



+ ^ [ Ki{\-sink)X%^{k)dk- ^ [ K, 



(41) 



+— / cosfc Xi(sinfc- A)X^ (A) dA. (42) 
Equations fiTj) and (|42|) lead to algebraic equations for the Fourier transforms of 6aN^{X) and 6pN^{k), which yields: 



^PNaik) = PNaik) — Poo Na{k) — -z— / ^1 (slu fc — A) COS (5(TAf^ ( A) rfA 

2n 



1. \ X- (k) — 

_^2u cosh [7r(A - sin k) /2u] ^" ^ ' 2tt 



dp, 



cosu;(A-/x) duj 

TT — p^exp(-wu)Xi (u) — . 

cosh(wu) ' ^■'^^'2n 



SpNa{k) = COSfc 

COS k 



doj 



+ 



2ti cosh(a;u) 

TT 1 



2n 



2u cosh [7r(sin k — X) /2u\ 



cos[Li;(sin k ~ sin g)]X^^ (g) dq 
X^JA) d\. 



(43) 



(44) 



Equations for 8u]ss^-i(X) and SpN^-2ik) are obtained in a similar way taking into account that in the thermodynamic 
limit kh = ±TT and A^ = ±00. This gives Eqs. (gT]), (gH) and (gS]), (g!]) where A^q is replaced hy Na - 2 and X^''^ by 
X'^''^. The quantity _2 is given by Eq. (|55|) in which the values of Xh are put equal to +00 and —00: 



Na/2-2 

^^„-2(A)-— E ^(A-A„)-aA,_2(A), 

^ a— 1 

and _2 by Eq. (j34p where the values kh are put equal to +7r and — tt: 



(45) 



J]] S{k - kj) + S{k -tt) + S{k + tt) 



PNa-2{k). 



The ground state energies (g} for the considered states can be rewritten as 
En, 



Na 

ENa-2 

Na 



= -2 



PNa{k) COS k dk — 2 / X 'f^ ^ (k) coa k dk, 

J — TT 

PN^-2{k) COS k dk — 2 / X^^_2(A:) cos fc dfc + — ^ cos fcft 



—2 / /9jVa-2(fc) cos fc dfc — 2 / X'ff^_2{k) cos k dk . 



The finite size corrections to the ground state energy are given by 

SEjM = En — Eoo,N, 



(46) 



(47) 



(48) 



(49) 
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where the energies Eao,N in the thermodynamic Umit are 

R — ly 

— — —21 Poo, Na{k) cos k dk, (50) 



r 4 

-2/ PooMa~2{k) COS k dk - 



(51) 



with Poo,w„(fc) given by Eq. (|35p and Poo,Na-2{k) by Eq. ([55)1 with fc/i = ±7r and A;i = ±oo. 

Using Eqs. (gZl), dH]), (46), dM]), and Eq. ^ with Poo,jv„ and /9oo,JV„-2 foUowing from Eqs. dSS]) and Eq. we 
obtain 



Na 

Na 



e,{k)X%(k) dk - e,{\)X%{\)d\ (52) 



e,{k)Xj,^^_^{k)dk- e,(A)X^_2(A) dA, (53) 



where 



e,(fc) = 2cosfc + 2 /°°:Mil£^I^L4!^cos(^sinfc) do., (54) 
Jq cj cosh(a;u) 

eJA) = 2 / , cos(ujX) dui, (55) 

Jq cosh(ijJM) 



and we used the fact that es(±cx3) = 0, so that Xf^ ^^^^ be replaced by Xf^ -2W Eq. ([55]) . 

For the repulsive Hubbard model at half filling, the gap is in the charge sector and the spin sector is gapless. 
Accordingly, the first term in the right-hand side of Eq. ([52|) describes the contribution of gapped charge excitations, 
and the second term is due to the contribution of gapless spin excitations. 

We now return to Eq. and using Eqs. ([55)1 write the finite size corrections to the gap of the attractive 

case in the form 

SA = SENa-2-SENa=SAng + SAg+S, (56) 

where 



SAng/Na^ I esi\)X'^^d\- I e,(A)X^^_2 rfA, (57) 
is the contribution of the gapless sector, 

6Ag/Na= I ec{k)X%^dk- I ec{k)X%^_^ dk, (58) 



is due to gapped excitations, and we used a relation 

2 



h—1 

The term S which is present in the gapped sector is given by 

d = J2e{k„)-e{7r)-e{-7T) (59) 

h 

and for Na Aoc ^ 1 it reduces to 

S^4{7r){k+-nr, 
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where k'^ is the value of kh which is close to tt at large Na- This term behaves as for large Na- In the limit 

of u <C 1 the value of k^ can be found from Eqs. ([57)1 and (1^0]) using the condition Zj^^_2 = 1/2 — 1/2 A'^^. For 
NaAoc ^ 1 we then obtain 

-^^n^] > (60) 



where at half filling the gap of the thermodynamic limit Aoo is given by 



wcoshtju 







at any interaction strength. As we will see below, the term S (j60p can become comparable with the (conformal) 1/Na 
correction for NaAoo ^ 40 and exceeds higher order power law corrections up to much larger values of Na- 

Note that the contribution of the gapped charge excitations for the repulsive Hubbard model corresponds to the 
contribution of gapped spin excitations for the attractive model, and the contribution of gapless spin excitations 
corresponds to the contribution of gapless charge excitations in the attractive case. 

The origin of the term 6 is the following. The state with two additional or two missing spin-up (or spin-down) 
particles of the initial attractive Hubbard model contains unpaired fermions with energies above the gap. Gapped 
excitations of the model are S = 1/2 - solitons which appear only in pairs. We thus calculate the exact energy of the 
states with 2 solitons which have different nonzero momenta and energies near the bottom of the excitation band. The 
contribution ((59)) takes into account these nonzero kinetic energies of the solitons and is proportional to the curvature 
of the excitation spectrum. 

It is worth noting that the term S is independent of the definition of the gap, except for minor changes for a 
relatively small number of particles. For example, it remains the same (as well as all other results of our work) if we 
consider the gap in the spectrum of triplet excitations: 2A E{N/2 + 1, N/2 - 1, -|[/|) - E{{N/2, N/2,-\U\). In 
this case we just get two unpaired spin-up fermions. 

The calculation for excited states of the system with unpaired fermions is the same as above. It simply assumes that 
kh are no longer fixed by the condition Zc{kh) = (Na — l)/2Na, but are related to the momenta ±p of the solitons. 
For the energy of the state containing two solitons with momenta p and —p we have e = 2A -I- 2[ec(7r — p) — ec(7r)], 
and for low- lying excitations at sufficiently large number of particles it is reduced to e = 2 A + e"(7r)p^, where 
A = Aoo + <5A„g + 6g. The minimum value of e is achieved at the minimum possible value of p which is equal 
to (tt — /c^), and we arrive at the gap A = A^o -I- SA, with 5A given by Eq. (|56|) and thus including 5 (|59p . It is 
important that this A is just the gap that is measured by the radiospectroscopy method used for obtaining the gap 
in experiments with two-component 3D Fermi gases Q. 



IV. EXPONENTIAL CORRECTIONS FOR HALF FILLING 



Throughout the paper we discuss the case where the inequality 



NaU > 1 



(62) 



is satisfied (except for Section VII), so that all analytical formulas obtained below remain valid for arbitrarily large 
Na- It is convinient to present the results in terms of the parameters NaAoo{u) and u because in the limit of u ^ 1 
the most important contributions to both SAg/Aoa and (5A„g/Aoo depend only on NaAoo (see below). 

In this Section we calculate the finite size corrections to the gap, originating from the gapped sector and given by 
Eq. (|551) . Using Eqs. p^ . (^5]) . (^5]) and the Poisson relation J2m ^(^ ~ '^") ~ X^m exppTrxmi] we obtain 



,{k)X$f{k) dk 



PNik) 



PN{k) 



cyi^[-2T:iNaZ%{k + iO)] + 1 eyi^[2T:iNaZ%{k - iO)] + 1 



dk. (63) 



To lowest order, we take the functions pN{k) and Z'ff{k) equal to their values in the thermodynamic limit. So, they 
are given by Eq. ^ and Eq. ((371) for the case N = Na, and by Eq. ^ and Eq. ^ for iV = TVq - 2, with kh = ±tt 
and Xh = ±oo. One can show that the integral from — tt to tt in the right -hand side of Eq. ([55)1 is equal to the 
integral from — tt + iarcsinhu to tt -I- zarcsinhu for the first term of the integrand plus the integral from— tt — zarcsinhu 
to TT — iarcsinhu for the second term. 
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FIG. 1: The correction 5Ag/Aoo versus NaAoo- The solid curve is the result of exact calculations from Eq. (|58p . and the 
dashed curve represents the result of Eq. (|64|l . In a) u = 1, in b) w = 0.5, and in c) u = 0.25 with the dotted curve showing 
the result of Eq. (|73p . with the replacement Na — > A^£i(l — u/tt) (see text after Eq. (|73|l ). The inset in c) shows the relative 
difference between the exact SAg and the exponential result of equation ([73} , k = (5A°'"'''* - SAl^^)/SA'g^'^^. 



Equation ([58)1 then reduces to 



Na 



PM^-2{k)ec{k) 



1 + exp [-2mNaZ'j^^{k)\ 1 + exp [-2TTiNaZ%^_^{k)\ 



dk 



1 



2mN„ 



4(fc)ln 



I + e^]y[2mNaZ%^{k)] 
1 + exppTTiTVaZ^ _2(fe)] 



dk + c.c. 



(64) 



where 7 — arcsinhu. 

It is convenient to present the ratio 5Ag/Ao 



as a function of A'^Aoo and u. In Fig. [T]we show (5Ag/Aoo versus 



A^aAoo for several values of u, and one clearly sees that for not very large this correction becomes significant. 

The edge points of the integration fco = ±7r ± iarcsinhtt are the saddle points at which pN^ (k) = dZ^^/dk = 0. For 
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sufficiently large Na we may use the saddle point approximation, and the expression for SAg becomes 

l4(fco)| 



dAa « C- 



TT^Na\p'{ko)\ 



exp[-5o] 



where 

4(fco) 

p'(fco) 
So 

c 

and we used the relation 



/■oo 

2u — 2\/ + 1 / Ji(a>) tanh(cjM) exp(— wit) 

"'0 



do; 



2^ 



- + {u^ + 1) / a>tanh(ciju) Jo{uj) exp(— wu) dcj 
Vu^ + 1 "'0 



7 ■ 



tanh(a;it) exp(— cjw) 



Jo{Lu)duj 



tanhfww) , , , \ 2i 
exp(— ww) aw = Zj^ (fco j H z-r- In 



The saddle point approximation assumes that the exponent in Eq. ()65p is large: 

So » 1. 

In the case of strong interaction, m ^ 1, Eq. (I65|) gives 



nN„, 



2r(5/4) 



r(3/4) 



Na-l ■ 



and one sees that in this limit the correction SAg is negligible. 
The situation changes for m < 1. In the limit of m ^ 1, from Eq. 



SA„ 



^j'^ ^oo exp[-A^A^j4] _ ^ 



we obtain: 
Aooexp[-AooA/'a/4] 



(65) 

(66) 
(67) 
(68) 

(69) 

(70) 
(71) 
(72) 



(73) 



where the gap in the thermodynamic hmit, Aoo, is given by Eq. pT|) . The criterion ([7T|) then becomes NaA^o ^ 1- 
The obtained relation (|73p is in accordance with the universal scaling behavior of the gap in massive quantum field 
theories [l^, which is expected for the Hubbard model at u ^ 1. For not very small u we should take into account 
corrections which are linear in u in the expression for So and in the preexponcntial factor for Aoo • This proves to be 
equivalent to the replacement Na Na(l — u/tt) in Eq. (|73p . Higher order corrections become important only for 
A^qAoo which are so large that the exponential contribution SAg is no longer important. 



POWER LAW CORRECTIONS 



The correction to the gap provided by the gapless sector, <5A„g, we calculate using the conformal field theory 1/Na 



expansion for the energy 



:ap proviaed 



En^ -Eo^ '^Vs{A+ + A_), Eo = eNa - -£^Vs, 



where eNa is the energy in the thermodynamic hmit, and the velocity of spin excitations is given by 

_ de e'^{B) 
' dp 2TTa{B) ' 

where B is defined by Eq. The conformal dimensions of primary operators are equal to 



(74) 



(75) 



(76) 
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80 100 



NaAoo 



FIG. 2: The correction (5Ang/Aoo versus NaA^o for u = 1. The sohd curve shows the resuh of exact calculation from Eq. (I57| 
and the dashed curve the result of Eq. ([83}. The dotted curve represents higher order power law corrections SAng (see text). 



Integer or half-odd integer numbers AiV^, Ds are characterizing the excitation states. The component of the charge 
dressed matrix S.s is determined from the equation 



6(A) = 1-^ J ^K2iX-v)Uv)dv- 



(77) 



For the zero-field case {B = oo) the solution is = 1/V2- The energy and momentum of a spin excitation are given 
by 



e.(A) = 2 



°° coswA ^ , .dto 

— Ji[uj) — , 

cosh WW oj 

°° sinwA .duj 
Jo{lj) — , 



cosh 6 



— ^ -2naiX). 



(78) 
(79) 
(80) 



The excitation state with the energy E{Na/2, Na/2 — 2) is characterized by the numbers Dg = 0, AA^^ = 1 
According to Eq. (fT4|) . this leads to the correction 



For the half-filled case we obtain 



_ 2^ /i(7r/2u) 
''^"^-7V,/o(7r/2u)- 

A more accurate result following from the calculations in Refs. [Ill, [l^ reads: 

_27r /i(V2w) / 1 

NaIo{7T/2u)\ 2\n[NaIo{TT/2u)] 



(81) 



(82) 



(83) 



The comparison of Eq. ((83)) with the result of exact calculations from Eq. ((57)l shows the validity of Eq. ((83)) even 
for not very large Na (see Fig. [2]). For example, at u = 1 even for Na = 10 {NaAoa ~ 13 ) the relative difference is 
~ 20%. The substraction of A„g ((83)) from the exact result of Eq. (I57p gives higher order power law corrections which 
we denote as SAng. For u—1 they are represented by the dotted curve in Fig.O 
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In the limit of strong coupling Eq. (|82p yields: 



^A„.-^; u»l, (84) 



and the power law correction ([5^ always dominates over the negligible exponential correction ([7^ . The situation 
changes for u < 1. For not very large NaAoc, the exponential correction SAg originating from gapped excitations 
becomes comparable with (5A„g and exceeds higher order power law corrections. We provide a detailed comparison 
of Ag with power law corrections in the next section. 




FIG. 3: The result of exact calculations from Eqs. (|57|l . (|58p and (|59() . In the left part the solid curve shows (5A„t,/Aoo versus 
A'^aAoo, the dashed curve SAg/Aoc, and the dotted curve S/Aoo- In the right part the solid curve shows the ratio of the higher 
order power law corrections (5A„g to the exponential correction SAg. In a) u = 1, in b) tt = 0.5, and in c) u = 0.25. 
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VI. DISCUSSION OF THE RESULTS 



We start the comparson of the non-conformal exponential correction to the gap with power law corrections in the 
limiting case of u ^ 1 and NaA^ 3> 1, so that Eq. ([75)1 is apphcable. Then, comparing the result of Eq. (175)1 
with that of Eq. (|55)l we see that the power law 1/Na correction dominates over the exponential correction for 
any A^qAqo significantly larger than unity. However, the exponential corection (|73)l is still larger than higher order 
power law corections in a wide range of Na- In the limit of Na — > oo the higher order corrections contain terms 
\ii[lii{Na)]/iNaln''{Na)), 1 / {Naln' (Na) , l/N^, etc., where g > 2 is an integer [ll|. Thus, for reasonable values of Na 
the term that should be compared with SAg ()73p is ~ ^/[Na In^(A^a)]- The argument of the logarithm may contain an 
w-dependent multiple B ^ exp(7r/2u) like the logarithm in the second term of Eq. ()83p . This is however not important. 
According to Eq. (fTT)l . in the limit of m ^ 1 the gap is exponentially small, Aoo ~ exp(— 7r/2u) <^ 1. Hence, for 
Na at which the exponential corection can still be important we have A^qAqo ^ and In(A^a) « ln(A;^^) ^ 1/u. 
So, irrespective of the presence of the w-dependent multiple B, the higher order power law corection becomes ^ 
(w^/iVa) ln(l/u) and it is smaller than the exponential correction ([73p for Na satisfying the inequality 

NaAoc<8\nQ-y, u<l. (85) 

Thus, at u <C 1 for NaA^o significantly larger than unity but still satisfying the condition ([85p. the exponential 
correction ([75)1 is legitimate and can be kept together with the power law correction (|83p. 

The situation is similar for intermediate values of u smaller than unity. This is seen from Fig. 3 where we present 
our numerical results or w = 0.25 and u — 0.5. However, already for u — 1 the exponential correction is smaller than 
the higher order power law terms (see Fig. 3) and should be omitted. As far as the term S is concerned, for u <C 1 it 
is comparable with the power law correction up to NaAoo ^ 40 and exceeds higher order power law corrections for 
much larger Na- Even if we compare S with the second term of Eq. ()83p . the latter is smaller at Na Aoo ^ l/w. 

In Fig. ID we present the results of exact calculations for the gap A at half filling from Eqs. ([2]), ([S]), (HI) and ITi)) for 
the same u as in Fig. |51 A direct comparison of A with the gap of the thermodynamic linit Aoo (|6ip shows that finite 
size corrections can be safely omitted for NaA^o > 40. On the other hand, already at NaAoc < 10, they dominate 
the gap. For u — 0.25 where Aoo = 5 • 10~^ this occurs already at Na < 1.5 • 10'^. For u = 1 we have Aoo = 1-28 and 
finite size corrections are important only for Na < 30. 

It is important that for NaA^o ~ 10 or even somewhat larger the non-conformal corrections originating from the 
gapped sector become comparable with power law corrections coming from the gapless sector. This is seen in Fig. |3| 
The exponential correction SAg is about 20% of (5A„g or smaller, and the correction d approaches (5A„g and can even 
exceed it. 

Qualitatively, the dependence A{Na) remains the same for smaller filling factors. This is seen from Fig. |5l where 
we present our numerical results for n = 0.2. For u = 1 finite size effects become important only at a very small 
number of lattice sites Na < 30. For u = 0.25 the thermodynamic-limit gap is Aoo ~ 6 • 10"^ and finite size effects 
are already important for Na « 500. 



VII. LIMIT OF Nau <C 1 



In the limit of NaU <C 1, which can be realized for u <C 1, the energy spectrum of the attractive Hubbard Model 
shows no exponential gap and both charge and spin sectors are conformal. The analysis of Lieb-Wu equations for 
this case has been done in [20l[2H. and found corrections ^ u/Na- This can be understood from the conformal 1/Na 
expansion as a consequence of the linear dependence of the velocities of elementary excitations on the interaction 
constant u 20, 22,]. 

Here we consider the limit of uNa ^ 1 for completeness and present first order corrections in u to the ground 
state energy and to the gap in the excitation spectrum. As in the previous sections, we calculate the energy for the 
Hubbard model in the repulsive case, where the Bethe Ansatz equations are easily solved, and then restore the energy 
for the attractive case using the particle-hole symmetry. 

So, consider a system of N particles {N^ spin-down and N — Ni spin-up) with repulsive interaction. From the 
Lieb-Wu equation ^ we obtain the momenta kj to first order in u: 

AT; Ni 

expUkiNa) = -^—r ^ => 5k, = fc,- - fc" = — > . — , (86) 

^ ' J-J- sin/c, - A„ - «M ■' ^ ^ iVa sin(/cO) - AO ^ ' 

a— 1 a— 1 ^ J ' " 
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FIG. 4: The gap A in units of t versus Na at half fiUing for it = 1 in a), = 0.5 in b), and u — 0.25 in c). The dotted line is 
the value of the gap in the thermodynamic limit. The sum 5Ang + SAg + 5 + Aoo (solid curve) coincides with the value of the 
gap found directly (crosses) from the Bethe Ansatz equations ([2]) and ((3| using Eqs. (jlj and (|14p . 



where fc" and A*^ are the momenta and rapidities for u —^ 0. The energy itself and the interaction-induced change of 
the energy are given by 



N 



■5-^ , n 4u v-A sin fc° 



N Ni 

A^a|J§sin(fcO)-AO' 



(87) 



where Eq is the ground state energy for u — > 0. We now calculate the densities of momenta k and rapidities A in the 
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Na 

FIG. 5: The gap A in units of t versus Na, calculated numerically for u = 0.25 in a), and u — 1.25 in b) for the filling factor 
n = 0.2. 



thermodynamic limit from Eqs. ([7]) and ([5]). To the lowest order in u we obtain: 

rQ 



2ct(A) = / (5(A-sinfc)p(fc)dfc, 

J-Q 
1 

p{k) = hcosfc / S{X — smk)a{X)dX. 

27r J_B 

The solution of Eqs. (HH) and dH]) is (nx < n^) 



Tini < k < 7r(n — n;). 



a(A) 



1 1 



2^ VT^A2 



(89) 

(90) 
(91) 



Then, using Eq. ([9]) for the total number of particles and the number of spin-down particles, we find an expression 
for the integration limits Q and B: 



B = sin (7rri|) , Q = 7r(n — rij). 
Using Eqs. (^0]) and (PT|) we obtain the interaction-induced change of the energy to first order in u: 

^ = 4u r dk r .'"'^ ■ ^(A)p(fc) dA = 



(92) 



Q j^gsink-X 
4w 4u r*"""^^ „ /"""^ sin^fc 

2 ' 



TT 



dk 







sin /c — sin q 



dq, 



(93) 
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where we turned to the variable q = arcsinA. For the case of attraction we should substitute n| ^ 1 — n^, n| ^ n^, 
in accordance with the symmetry properties (|12p . Integrating over dq we find: 

5E{-u) ^ . 2 2u / , , /tanfc + tanTrn, \ „ 

^ ' - -4wix + ^unl + ^ 2 / - / tanfc In ^ dfc, (94) 



Na ^ 7r2 I / \tan/s - tauTrnx 



where n| and are already occupation numbers for the attractive Hubbard model, and we assume that < n|. 
Eq. ([94]) leads to the following result for the interaction-induced change of the energy to first order in Nau: 



SE(—u) , 4u u , / 1 + tan^ TTfit \ , /tanTrnt + tauTrn 



= 4ur7,i — 4un| H arctanvrnj 2 -'^'^ 2 ) 



Na ^ TT TT^ \l + tan Trnj/ \tan7rnf — tauTrnj^ 

2u / „ / 2tan7rni \ „ /tanTrnt + tanTrni \ „ / tan Trot — tan Trrt 
' 'ReLi2 — - ReLi2 — - + ReLi2 ' ' 



TT^ \ \tan7rn|— \ tanTmj^— i J \ —tamrni—i 

where Li2 (z) = X^fcLi /k"^ is a polylogarithmic function. 

In the limit of small filling factors, ^ Na and <^ Na, after a straightforward algebra we obtain: 



SE 



u(-4nx - 4n-rnx + 2n^) (95) 



The limit of small filling factors in the Hubbard model corresponds to the gas phase of spin-1/2 fermions. For this 
case the ground state energy at Nu <C 1 has been calculated in Refs. [H, [U13, and the result of Eq. (|95p coincides 
with that of Refs. in the attractive case. 

Using Eq. ([6]) we then find a small interaction-induced correction to the gap in the excitation spectrum (n| = — 
n/2) of the attractive model to the lowest order in uNa- For small filling factors we have: « Au/Na, and in the 
considered limit of NaU ^ 1 this correction is small compared to the level spacing ~ 1/Na in our finite size system. 



VIII. CONCLUSIONS 



In conclusion, we have studied finite size effects for the gap in the excitation spectrum of the ID Fermi Hubbard 
model with one-site attraction. For the situation in which the thermodynamic-limit gap A^o exceeds the level spacing 
(near the Fermi energy) of the finite size system, there are two types of finite size corrections. For large interactions 
(m S> 1) the leading is a power law conformal correction to Aqo, which behaves as 1/Na and originates from the gapless 
sector of the excitation spectrum. We also find non-conformal corrections originating from the gapped branch of the 
spectrum. As found at half filling, in the weakly interacting regime (u < 1) the non-conformal corrections can become 
of the order of the conformal correction even for the number of particles (lattice sites) as large as ~ 20/Aoo. Also, for 
M <C 1 and large Na Aoo, the exponential correction ([75]) is legitimate as long as the condition ([55]) is satisfied. Thus, 
we have the full right to take it into acount together with the power law corection ([55]) . 

For sufficiently small number of lattice sites (particles) the gap A is dominated by finite size effects. From a general 
point of view, this happens when A^c £ 1/Na, i- e. Aqo is smaller than the level spacing of the finite size system at 
energies close to the Fermi energy. Accordingly, for large interactions (m 3> 1) the finite size effects are not important 
as long as Na ^ 1. However, in the weakly interacting regime (u < 1) they become dominant already at significantly 
larger Na than a simple dimensional estimate I/Aqq. This is clearly seen from our results in Fig. 3 and Fig. 4 for 
A{Na) at haff fiUing. 

Our findings are especially important for the studies of the ID regime with cold atoms, where the number of 
particles in a ID tube ranges from several tens to several hundreds [1, For such a system in the weakly interacting 
regime one can not use the result of the thermodynamic limit for the gap. Consequently, one can not employ the 
local density approximation for A based on this result, for finding the spectrum of isospin gapped excitations in an 
external harmonic potential. 
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